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We present a field-theoretic renormalization group analysis of a polymer chain 
immersed in a binary good solvent close to its critical demixing point. We first show 
that this problem can be mapped on a bicritical field theory, i.e. a ($ 2 ) 2 -model with 
a mass anisotropy. This implies that the end-to-end distance of the polymer is now 
controlled by a new critical exponent vb related to the quadratic mass anisotropy 
operator B. To show this we solve the RG equation and calculate explicitly the 
exponents and the mean end-to-end length of the chain. 
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In this paper we consider the critical behavior of a polymer chain in a critical fluid pL This 
problem is quite poorly understood, because its a coupled system with two macroscopic length 
scales. The physics is determined by an interplay of the chain length and the correlation length of 
the critical fluid, which may undergo a phase separation at a certain temperature. The polymer in 
solution and the critical fluid are separately successfully described by scaling and renormalization 
group theories. The separation of microscopic and macroscopic length scales leads there to uni- 
versality. That means the main properties of the system are independent of microscopic details, 
which only change prefactors in the physical quantities and make analysis of experimental results, 
at least in their limiting behavior straightforward. We will show that also for our coupled system 
the application of the RG is possible and yields a description of the conformation of the chain. 
First insight beyond mean field theory was brought in by Brochard and de Gennes by the use 
of scaling theory and by the simulation of Magda et al ||, but still there is a lack of a quantitative 
investigation of the chain conformation near the consolute point of the solvent. The start point of 
such an investigation is the settled up field theory by Vilgis et al |Q . 

The polymer chain of configuration R(s) is modeled by an Edwards- Wiener functional [p| 



where s is the dimensionless contour variable, I the Kuhn segment length, d the space dimension, v 
the parameter of the bare short ranged excluded volume interaction, TV is the number of segments 
and L = IN the contour length of the chain. The reduced Hamiltonian of the fluid can be given in 
terms of three short range pair-potentials between the good solvent molecules 

<T,T i,j 

where a, r = 1,2 and i,j label the molecules. The interaction between the chain and the fluid is 
described as 



Hi = [ dsJ2J2 v m,S(R(s) - rj) , 

J ° a=1.2 i 



(3) 



where v ma characterizes the short range excluded volume interaction between the monomers and 
the different species of the fluid. If their solvent quality is slightly different, that means v m \ > v m i 
or v m \ < v m 2, there is preferential affinity, the chain is likely to be surrounded by solvent 2 or 
solvent 1. After the use of collective density fields |I| for the fluid and the chain we get the equivalent 
complete field theoretical Hamiltonian 
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n = Jd d r |i(V0(r)) 2 + \y^ 2 (r) + \(^(r)) 2 + l -r ^{r) + 

+ ^ 2 (^ + |(0 2 (r)) 2 + |v/M+2|0 2 (r)V, 2 (r) . (4) 

The polymer field <j> is an M component field with M — » after the evaluation of the perturbation 
expansion. The quantity 4> 2 is related to the monomer density. This follows from the well known 
de Gennes transformation. Whereas ip is a scalar field which describes the relative concentration 
fluctuations of one of the fluid constituents. The term <f> 2 tfj 2 originates from the incompressibility 
condition p| and is important for the physics: In the case of preferential affinity (wq > 0) it re- 
swells the chain again in the critical region of the consolute point ^,|| . Whereas in the case of no 
preferential affinity it is the lowest order coupling between fluid and polymer. It is interesting to 
note that the same Hamiltonian but with an imaginary coupling wq appears in the modeling of a 
single screened polyelectrolyte chain, i.e. a Debye Hiickel chain with corrections ||. 

In the following we consider the case with no preferential affinity, i.e. the cubic coupling wq = 0. 
We may write the Hamiltonian as one with a composite M = M + 1 component field, with M — ► 1 
due the polymer limes. The fixed point of such a system is for M < 4 the Heisenberg fixed point 
with all corresponding renormalized quartic couplings equal u* — u\ — u\ — u*{M) 0. Due 
to universality we are allowed to set these couplings already equal in the original (renormalized) 
Hamiltonian |?J and we get 

n = J d d r |i(V*(r)) 2 + \y^ 2 {r) + \t ^{t) + ^(* 2 (r)) 2 | , (5) 

with <& = (<fi, ip). This is exactly the Hamiltonian of a bicritical field theory and well investigated by 
Amit and Goldschmidt ||, but now with yo as a monomer chemical potential, an inverse Laplace 
variable. To analyze the general behavior of the end to end distance, we use the so called soft 
expansion. That means we write (|^) in the form 



H = J d d r|i(V*(r)) 5 



l -t^(r) + ^(r)f-^B(r)) . ((,) 



with the definitions 



B(r) = 1/M ((M - l)^ 2 (r) - 2 (r)) 
t Q = 1/M(r + (Af - l)y ) 

9o = t - y . (7) 

The advantage of this splitting is that <fr 2 and B are orthogonal under renormalization, that 
means divergences which are generated by $ 2 (B) insertions in vertex functions are absorbed 
completely by the corresponding counterterms of 3> 2 (B) and these counterterms do not mix 
with each other under renormalization. We calculate the new exponent vb which is given by the 
anomalous dimension of the operator B corresponding due to the limit M — > 1 to the monomer 
density. The other exponents rj, v and the /3-function are the usual Heisenberg ones, which is 
obvious in dimensional regularization Q (the masses go and to are then allowed to set to zero). 
The upper critical dimension d c — 4 follows from a dimensional analysis of as for the symmetric 
Heisenberg model with go — 0. . We show that 

v B = 1/(2 -7b) 

= (l + e /9 + 0(6 2 ))/2 , (8) 

and the cross over exponent is 

4> = v/v B 

= l + e/18 + 0(e 2 ) , (9) 
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and v = (l + e/6 + 0(e 2 )/2, consistent with ||. But the soft expansion (||) is only useful for general 
considerations like the solution of RG-equations. Whereas for perturbation expansions of polymer 
quantities it is not possible to use (^), because an expansion in is an expansion in an inverse 
Laplacian and therefore all back transformed contributions are order by order zero. 

Now we perform a general renormalization group analysis of our model given by equation (^|). 
Our goal is to calculate the end to end distance of the polymer 

<R 2 > = -T c 1 (d q .& 2 \ q )\ q=0 )) /TJ 1 (&*>(q = 0)) 

= D(u)F(r,C) , (10) 

here is G^(r) = J d d q/ {2n) d e l * r &V ( q ) = M J T><$>Vij} §(r)$>{Q)e- H ^\ the propagator of 
the polymer, Af is a normalization constant, is the inverse Laplacian with respect to (wrt) 
C = 21L = 2l 2 N/d, & 2 '(q) is the Fourier transform of the renormalized chain propagator and 
D(u) is a non universal constant. 

Due to the particular choice of B we have a separate (diagonal) renormalization scheme 0, 
which expresses the bare parameters in terms of their renormalized ones: 

t a = Z$st, g = Z B g, u Q = Z u uSd^ (11) 

Here /_t _1 is an external length scale and Sd = VLd/{2ir) d with the volume of the unit sphere. 

We denote an N point vertex function as T^ N ' and a composite one with N points and L, $ 2 and 
K, B insertions as r (N < L > K \ The wave-function renormalization factor is calculated by expressing 
the renormalized vertex functions in terms of the bare ones as 

rf'^fr^' . (12) 

Whereas the renormalization of the composite ones give similar Q 

T (N,L,K) =Z m Z L aZ K T (N t L,K) (13) 

The renormalization group equations follow from an variation wrt the external scale \x and are 
given by 

+ I3d u + 7*2 rd T + lB gd g - 7 *iV/2] rjf' = , (14) 

with the Wilson- functions (3 — ^idu/dfi, 73 /4> 2 = ~P(d In Zg /$2 / du) and 7$ = (3{d In Z$ / du). The 
equivalent form in polymer variables is obtained by an inverse Laplace Transformation with partial 
integration wrt g. The vertex functions are calculated using a dimensional regularized scheme with 
minimal subtraction of poles 0. The vertex function yields the wave function renormalization 
Z<p = 1 + 0(u 2 ). The determination of Z u is obtained by the 1/e-part of T^ 4 ' and is 

37/ 

Z u = l + ^ + 0(u 2 ) , (15) 

from which follows the /3-function the zero of which yields the infra-red stable fixed point u* = 
2e/3 + 0(e 2 ). An evaluation of F^ 1 ' ) yields 



whereas F^ 2,0 ' 1 ) yields 



Z^=l + - + 0(u 2 ) , (16) 



Z B = l + - + 0(u 2 ) . (17) 



From the above equations we now obtain the exponents due to the relations 



3 



?7 = 7*OV), v 1 - 2 = 7$ 2 (u*,e), i^ 1 - 2 = 7b(w*, e) , (18) 

which yield the exponents given in equations (||) and (|^). Our results are to order e, with corrections 
of the order of e 2 . 

To see the scaling form of < R 2 > we need to solve the RG equation for the Fourier transformed 
chain partition function Z R 2 \q) = T^ 1 (& 2 ^(q)) which is given by 

[fid/j, + pd a + 7$2T<9 T - jbL8 l +7$ +75] Z R ' (q,r,L,u, n) = , (19) 

The solution is obtained by the method of characteristics fjj and gives the scaling form of < R 2 > 
to be § 

F(t,C) = C 2ub F{tC ub/v ,u*) 

= \t\- 2u F ± (jC\t\ u / Vb ,u*) (20) 

We see that in the above equations the cross over exponent <fi is appeared. Here ± denotes T > T c 
and T < T c , respectively. 

Now we calculate the explicit form of the scaling function of the end to end distance < R 2 > 
in case of r — > (T > T c ) and C — > 00 so that tC stays finite. We need the general matrix 
renormalization scheme, because no soft expansion wrt the monomer chemical potential is possible. 
We have to use instead of equation ([n]) for the vertex renormalization 

yo = z 2 it + z 22 y 

t = Z\\r + Z 12 y . (21) 
These Z-factors are simply given by the renormalization of the corresponding vertex functions with 

~ f2) ~ (2) 

zero momentum for the polymer and the fluid T 22 (k — 0) and (k = 0) which fix the values of 
Z21 and Z22 and Z\\ and Z\ 2 . They are given to order 1-Loop in the figures 1 and 2. 

1 

o 

r 22 (fc) = - 2 g + 2 ' ' 2 + 2 1 — 2 + °(" 2 ) 

FIG. 1. The diagrams for the polymer V 2 2 ■ A line with the number 2 cor- 
responds to a free polymer propagator (k 2 + y)^ 1 , a line with 1 to a fluid 
propagator (k 2 + t)" 1 and doted line to a interaction vertex u/4\. Note that 
we do not include external legs and evaluate the graphs at k = for the ver- 
tex function. The chain propagator G 1 - 2 ' is also calculated using these graphs 
including external legs and with k finite. 
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o 

fnW = 1 i + j— ^ — i + 1 : 1 + °(" 2 ) 

FIG. 2. The diagrams for the fluid T<$ 



The evaluation of their Feynman diagrams yield Z 2 i = £ + 0(u 2 ), Z 2 2 = 1 + 57 + 0(u 2 ), 
Z X x = 1 + Ye + 0{u 2 ) and Z 12 = 0, in the polymer limit M -> 1. 

~ (2) 

The corresponding RG equation for ^ 22R is then 

\fidn + (3d u + K T rd T + (k 2 it + K 2 2y)d y - 7$] T 22 ' R = , (22) 

and follows from a variation wrt the external length . The new coefficients are k t = —/3d u In Z\\, 
K22 = -/3d u lnZ 2 2 and k 2 i = -l/Z 2 2(0d u Z 2 i - Z 2 ik t ). 
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After an inverse Laplace transformation wrt y of the RG equation for the propagator = 
l/r^ 2 ) we get for the 2-point chain partition function 

[fidf, + /3d u + TK T d T - k 2 2L8l + 7* - ^22 - k 2 itL] Z 22 ' r = . (23) 

With the use of characteristics as mentioned after (|f9| ) we find the solution of (|23| ) and it turns 
out the same scaling form for < R 2 > as in the soft expansion (^) because both renormalization 
schemes are equivalent. But the critical exponents are now defined as \/vb = 2 — H22{u*), l/v = 
2 — k t (u*) and r\ = 7$ (it*). Because the solution of the general vertex matrix renormalization 
has to be consistent with the diagonal (soft expansion) scheme of (|l9|), it follows that Z22 = Zb 
and Zn = Z$,2 to all orders in u. But we need the matrix scheme for the analysis of the explicit 
loop-expansion: By using the flow equations from the method of characteristics similar as in the 
calculation of the equation of state for magnetic systems we get from the expansion of f ^ 
finally the explicit result valid to order e 

< R 2 > ~ L 2VB . (24) 

This scaling form is consistent with the general result (]2p|). 

Now we express our renormalized quantities by the original bare ones and find that 

< R 2 >= D{u)l 2 N 2uB , (25) 

with a non universal constant D(u) There is in 1-loop order no explicit temperature dependence 
of the demixing binary solvent but we have the implicit condition that tqNI 2 is finite in the critical 
limit. 

To summarize: we found that in the critical region the demixing tendency of the fluid yields, 
if the number of monomers N — > 00 and the critical temperature deviation of the fluid r — > 
such that INtq stays finite, to a smaller chain size than in each of the good solvent components 
separately. This is because the new exponent vb is less than the self avoiding walk exponent 
b>o = (1 + e/8 + 0(e 2 ))/2. To get an explicit demixing temperature dependence in the end-to-end 
distance of the polymer chain a calculation to a higher Loop order is necessary. Now only the 
exponent vb can be found in 2-Loop order in reference 0. 

The binary solvent with preferential affinity to the monomers is not tractable by the renor- 
malization group. Due to the induced preferential adsorption the chain likes to be surrounded by 
the better solvent. This effective long range attraction of the range of the fluid correlation length 
yields in a specific temperature range to a partly collapsed chain [|] ||| . In the RG calculation this 
manifests in a dominant cubic term in the Hamiltonian which implies an imaginary solution of the 
(3- function M. 

A related problem is a screened polyelectrolyte chain in a demixing binary solvent. The effective 
attractive interaction induced by the preferential affinity leads then to an effective charge of the 
monomers 0]. 

Real systems have preferential affinity. Therefore its interesting to investigate if at least right at 
the consolute point of the fluid the preferential adsorption term has still significant influence on 
the (infinite) chain size, because the chain is then again swollen. This could be done by computer- 
experiments. 

Note added: After acceptance of the present paper an experimental study appeared [ |To[ which 
agrees in several points with our results. 

We acknowledge helpful discussions with H. K. Janssen, T. M. Liverpool and J. J. Magda. We 
particularly thank H. K. Janssen and T. M. Liverpool for a critical reading of the manuscript. 
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o 

f *(*) = 5 1 + 2 ; ' 2 + i =—5 + 0(« 2 ) 

FIG. 3. The diagrams for the polymer r 22 . A line with the number 2 corresponds to a free polymer 
propagator (k 2 +y)~ 1 , a line with 1 to a fluid propagator (k 2 + t) _1 and doted line to a interaction vertex 
w/4!. Note that we do not include external legs and evaluate the graphs at k = for the vertex function. 
The chain propagator is also calculated using these graphs including external legs and with k finite. 
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rn(*o = 1 i + 1 — ^ — i + 1 — ; — i + °( u2 ) 



FIG. 4. The diagrams for the fluid rff. 
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